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I.   INTRODUCTION 

Although  it  is  well  recognized  that  combat  is  a  complex 
random  process,  a  model  that  provides  "closed-form" 
solutions  can  often  be  an  aid  to  the  military  analyst  by 
explictly  portraying  interrelationships  between  combat 
parameters.  Deterministic  differential  equations  are 
frequently  used  to  model  the  combat  process,  and  vary 
greatly  in  their  detail  and  level  of  complexity.  The 
military  analyst  must  be  able  to  recognize  significant 
variables  and  their  relationships.  An  idealized  model  is 
useful  for  giving  a  "first  cut"  portrayal  of  such 
relationships.  In  particular,  one  can  often  obtain  an 
analytical  solution  to  such  models.  An  analytical  solution 
that  can  be  expressed  in  terms  of  "elementary"  functions  can 
often  provide  insights  into  the  combat  process  and  identify 
factors  affecting  battle  outcome. 

The  use  of  differential  equations  to  model  warfare  was 
pioneered  by  F.  B.  Lanchester  in  1914  [Eef  8].  This  original 
work  was  quite  idealized  and  assumed  that  the  fire 
effectivenesses  were  constant  throughout  the  battle.  Many 
extensions  to  Lanchester1 s  classic  formulation  have 
subsequently  been  developed,  to  include  such  factors  as 
mobility  of  forces  and  fire  effectiveness  that  varies  during 
the  battle.  Taylor  and  Brown  [Ref  13]  have  developed  a 
mathematical  theory  for  solving  Lanchester-type  equations 
with  these  temporal  variations  in  fire  effectiveness.  The 
fire  effectiveness  of  forces  is  represented  by  what  are  now 
called  Lanchester  attrition- rate  coefficients.  Taylor  and 
Brown  give  results  for  two  specific  forms  of  the 
coefficients:  (1)  effectiveness  of  each  side's  fire 
proportional  to  a  power  of  time,  and  (2)  effectiveness  of 
each  side's  fire  linear  with  time  but  a  nonconstant  ratio  of 
attrition-  rate  coefficients. 


Taylor  and  Brown  present  accurate  numerical  solutions 
to  equations  with  the  above  attrition-rate  coefficients. 
These  solutions  are,  in  general,  stated  in  terms  of  infinite 
series.  The  solution  method,  however,  is  rather  complicated 
and  involved  the  development  of  new  mathematical  functions 
alien  to  the  military  analyst.  Exact  results  are  certainly 
beneficial,  but  the  analyst  also  needs  a  solution  in  terms 
of  "familiar"  functions  to  permit  parametric  analysis  and 
simplify  recognition  of  significant  variables.  To  provide 
such  a  solution  form,  Taylor  [Ref  14]  has  suggested  the  use 
of  the  so-called  Liouville-Green  approximation,  which  can  be 
expressed  in  terms  of  relatively  simple  functions.  For  an 
approximation  to  be  useful  in  predicting  force  levels, 
however,  it  must  be  reasonably  "close"  to  the  exact 
solution. 

This  paper  provides  a  numerical  evaluation  of  the 
Liouville-green  approximation  for  the  two  forms  of 
attrition-rate  coefficients  mentioned  above.  A  comparison  is 
made  between  the  approximation  in  each  of  these  cases  and 
the  numerical  results  given  by  Taylor  and  Brown .  Numerical 
examples  and  analytic  considerations  will  be  utilized  to 
investigate  the  accuracy  of  the  approximation  in  these 
cases. 

The  organization  of  this  thesis  is  as  follows.  We  first 
review  Lanchester's  classic  model  of  combat  between  two 
homogeneous  forces.  These  results  are  then  extended  to  the 
variable  attrition-rate  coefficient  case,  and  the  specific 
forms  to  be  studied  are  discussed.  The  representation  of 
the  solution  to  these  variable-coefficient  equations  in 
terms  of  so-called  general  Lanchester  functions  is 
presented.  Then  Liouville's  normal  form  and  approximation 
are  introduced,  evaluated,  and  discussed  for  the  two 
specific  forms  of  attrition-rate  coefficients. 


II.   LANCHESTER'S  CLASSIC  FORMPLATION 

F.  8.  Lanchester  was  an  English  aeronautical  engineer 
who  first  developed  a  mathematical  formulation  for  combat 
between  two  conflicting  forces.  He  published  this  theory  in 
a  series  of  articles  in  the  British  journal  Engineering  in 
1914  £Bef  8].  Lanchester1 s  objectives  were  to  provide 
insights  into  "modern"  combat  and  to  quantitatively  justify 
the  military  principle  of  concentration.  He  stressed  that 
advances  in  technology  had  brought  long-range  delivery 
capability  of  modern  weapons  which  allowed  for  concentration 
of  firepower.  Lanchester  hypothesized  that  under  these 
modern  conditions,  combat  between  two  military  forces  could 
be  modelled  by: 


dx/dt  =  -ay      with    x  (t=0)  =  x 

0 


/ 


dy/dt  =  -bx      with    y  (t=0)  =  y   ,        (1) 

where  x  (t)  and  y (t)  are  the  numbers  of  X  and  Y  forces  at 
time  t,  and  t=0  denotes  the  time  at  which  the  battle  begins. 
These  equations  are  valid  only  for  x  and  y  greater  than 
zero.  The  parameters  a  and  b  are  non-negative  constants  that 
have  come  to  be  referred  to  as  Lanchester  attrition-rate 
coefficients.  These  coefficients  represent  each  side's  fire 
effectiveness  or  firepower.  The  most  common  set  of 
circumstances  under  which  these  equations  have  been 
hypothesized  to  apply  is  that  both  sides  use  aimed  fire  and 
target  acquisition  times  are  constant  [Ref  16]. 

From  (1)  Lanchester  deduced  his  classic  square  law 

2    2  2     2 

b(XQ  -  x  (t))  -  a(y   -  y  (t) )  ,  (2) 

which  has  the  important  implication  that  a  force  can 
significantly  reduce  its  casualties  by  initially   committing 
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more  forces  to  the  battle.  The  square  law  also  shows  that  in 
this  simple  model  the  Y  force  wins  (ie.,  X  is  annihilated) 
if  and  only  if 


x  /y  <  Va/b"  .  ( 3 ) 

0   0 


From  (1)  Lanchester  also   derived  the   now   well-known 
results  for  the  time  history  of  the  force  levels 


x  (t)  =  x  coshv/ab  t  -  y  Va/b  sinh\/ab  t  , 
0  0 

y  (t)  =  y  cosh\/a~b  t  -  x  v/b/a  sinh\iaFt  ,      (4) 
0  0 


The  victory-prediction  condition  (3)  was  given  for  a 
"fight  to  the  finish".  As  H.  K.  Weiss  [ Ref  15]  points  out, 
most  battles  terminate  before  complete  annihilation  of  one 
force  occurs.  Different  models  of  battle  termination  have 
been  proposed,  among  them  that  combat  ends  when  either  of 
two  given  "breakpoint"  force  ratios  is  reached.  Taylor  has 
shown  that  (3)  is  necessary  and  sufficient  for  I  to  "win"  a 
fixed-force-ratio-breakpoint  battle  [Ref  14]. 

Lanchester's  original  formulation  (1)  involved  many 
implicit  assumptions  that  limit  its  usefulness  for  modelling 
real  combat.  The  forces  were  assumed  to  be  homogeneous, 
there  were  no  replacements  or  withdrawals,  and  no  movement 
of  forces  was  considered.  Setting  the  attrition-rate 
coefficients  a  and  b  in  (1)  as  constants  implies  that  the 
fire  effectiveness  per  firer  of  each  side  is  constant 
throughout  the  battle. 
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III.  VARIABLE  ATTRITION-RATE  COEFFICIENTS 

When  modelling  mobile  weapon  systems  and  movement  of 
forces,  one  recognizes  that  the  assumption  of  constant  fire 
effectiveness  is  open  to  question  [Ref  2].  If  fire 
effectiveness  is  permitted  to  vary  during  the  battle, 
Lanchester*s  formulation  (1)  takes  the  form: 

dx/dt  =  -a(t)  y     with  x  (t  =  0)  =  x   , 

0 

dy/dt  =  -b(t)  x     with  y  (t=0)  =  y   ,        (5) 

where  a  (t)  and  b (t)  are  time-dependent  attrition-rate 
coefficients.  The  form  of  the  attrition-rate  coefficients 
depends  on  such  variables  as  force  separation  (range  of 
battle) ,  tactical  posture  of  the  targets,  firing  rate  and 
rate  of  target  acquisition  [Ref  5]. 

New  operations  research  techniques  for  predicting  these 
attrition-rate  coefficients  were  developed  in  the  1960's.  A 
significant  contribution  in  this  area  was  the  development  of 
a  methodology  by  S.  Bonder  [Ref  2  and  4]  for  prediction  of 
coefficients  based  on  weapon  system  performance  data. 
Another  important  development  was  the  method  of  G.  Clark  for 
the  maximum  likelihood  estimation  of  these  coefficients 
based  on  aonte  Carlo  simulation  output  data  [Ref  6].  The 
work  of  these  and  others  has  facilitated  the  use  of  models 
such  as  (5)  in  current  combat  studies. 

To  effectively  model  today* s  highly  mobile  warfare,  the 
form  of  the  attrition-rate  coefficients  in  (5)  must 
accomodate  current  tactics  and  battle  conditions. 
Specifically,  the  coefficients  must  permit  the  fire 
effectiveness  of  one  or  both  sides  to  be  non-zero  at  the 
start  of  battle  (to  model,  for  example,  an  ambush).  The 
coefficients  must  also  allow  for  the  opposing  weapon  systems 
to   have   different   maximum  effective  ranges.   Since  mobile 
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forces  close  with  each  other  in  combat,  their  fire 
effectiveness  should  be  permitted  to  increase  or  decrease  as 
a  function  of  range  or  time.  The  attrition-rate  coefficients 
examined  in  this  paper  have  the  flexibility  to  model  these 
various  situations.  Motivated  by  previous  work  by  Bonder 
and  Farrell  [Ref  5],  Taylor  and  Brown  have  considered  the 
following  attrition-rate  coefficients: 

m  n 

a(t)  =  k  (t+C)     and   b  (t)  =  k  (t+C+A)  .     (6) 

a  b 

The  specific  cases  of  (6)  to  be  considered  here  will  be  that 
in  which  A=Q,  C>0, 

m  n 

a  (t)  =  k  (t+C)     and   b  (t)  =  k  (t+C)  ,      m 

a  b  v '  ' 

and  that  in  which  m=n=1,  A>0,  C>0, 

a(t)  ■  k  (t+C)     and   b  (t)  =  k  (t+C+A)  .      (8) 

a  b 

The  parameter  A  will  be  referred  to  as  the  "offset 
parameter"  since  it  permits  modelling  of  combat  between 
forces  whose  weapon  systems  have  different  maximum  effective 
ranges.  The  parameter  A  also  allows  the  ratio  of 
attrition-rate  coefficients  to  be  non-constant.  Koopraan  [see 
Ref  10]  first  observed  that  if  the  ratio  of  attrition-rate 
coefficients  is  constant,  the  time  solution  to  (5)  is  no 
more  complicated  than  the  constant-coefficient  solution. 
That  is,  when 

a(t)  =  k  h  (t)     and    b  (t)  =  k  h  (t)  ,       (9) 

a  b 

where  h (t)  is  the  common  time-dependent  factor  in  both 
attrition-rate  coefficients,  the  solution  is 


x  (t)  =  x   cosh  K  (t)  -  y  Jk  /k   sinh  w(t), 
0  0   a   b 
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where   W  (t)   =  \j  k  3c  J   h(s)  ds.    Koopman  also  found  that  the 

a  b  o 

constant  ratio  of  coefficients   again   yields  a   square-law 
relationship (see  also  Ref  16). 

The  parameter  C  in  (7)  and  (8)  allows  for  the  initial 
fire  effectiveness  of  each  side  to  be  non-zero.  This  permits 
modelling  of  battles  which  begin  within  the  maximum 
effective  ranges  of  the  opposing  weapon  systems.  The 
parameters  m  and  n  allow  the  fire  effectiveness  of  each  side 
to  increase  or  decrease  with  time.  With  m,n<0,  the  fire 
effectivenesses  decrease.  This  situation  might  model  a 
battle  in  which  combatants  take  cover  and  improve  their 
positions  as  the  battle  progresses.  With  mrn>0,  the 
effectivenesses  increase,  such  as  might  be  the  case  in  which 
two  forces  are  closing  with  each  other. 

To  model  an  engagement  in  which  forces  close  with  one 
another,  it  is  beneficial  to  relate  their  instantaneous  fire 
effectiveness  to  the  range  separating  them.  As  this  range 
decreases,  one  might  expect  fire  effectiveness  to  increase. 
The  attrition-rate  coefficients  (6)  can  be  expressed  in 
terms  of  range  between  forces  and  the  maximum  effective 
ranges  of  the  opposing  weapon  systems: 

m 
dx/dt  =  -a(r)  y  =  -a  (1-r/2  )  y, 

0      a 

dy/dt  =  -b(r)  x  =  -b  (1-r/R  )  x.  (10) 

0       b 

These  equations  have  been  used  by  S.  Bonder  [Refs  1   and   3] 

to   model   a   constant-speed   attack   on   a  static  defensive 

position.   Figures  1  and  2  show  sample  relationships  between 

the   attrition-rate   coefficients  and  range  between  opposing 

forces  for  various  input  parameters.   In  these  equations,   r 

is   the   range   between   the   forces,   and  R   and  R   are  the 

a      b 
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maximum  effective  ranges  of  the  1     and  X   weapon   systems, 

respectively.   The  constants  a  and  b   represent  the  fire 

0      0 

effectiveness  of  each  side   at   r=0.   Range  and   time  are 

related  by  r  (t)  =R  -vt,  where  R  denotes  the  opening  range  of 

0  0 

battle,    and   v      is      the     constant     attack      speed.      With      this 

relationship,   the   attrition-rate   coefficients  in    (10)    can  be 

related   to    (7)    and    (8) .    For    (7)    we  have: 

k     =  a    (v/R  )m,     k     -  B    (v/R  )n,  c  =    (R  -H)/ti(11) 

aOa  bOb  aO  \ill 

where  R   =  R   >  R  ,  while  for  (8)  we  have: 
a     b    0 

k      =    a    v/R      ,    k      =    B    v/R    ,       A    =    (R    -R    ) /v,      C   =     (R    -H    )/v, 
aOabOB  ba  aO 

where   R      >   R      >   R    .  ' 12^ 

b  a  0 
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IV,  A  STANDARD  FORM  FOR  THE  VARIABLE-COEFFICIENT  SOLUTION 

The   variable   attrition-rate  coefficient  equations  (5) 
yield  the  force  level  equations 

2    2 
d  x/dt  -  [1/a(t)  da/dt]  dx/dt  -  a(t)b(t)  x  =  0,  (13) 

with  initial  conditions 

x(t=0)  =  x    and   [  1/a  (t)  dx/dt  ]    -   -y  , 
0  t=0     0 

and  similarly 

2     2 

d  y/dt  -  [1/b(t)  db/dt]  dy/dt  -  a(t)b(t)  y  -  0r  (14) 

with  initial  conditions 

y  (t=0)  =  y    and    [  1/b  (t)  dy/dt]    «  -x  . 
0  t=0     0 

The  solution  to  (13)  is  given  by 

x(t)  =  C  X  (t)  +  C  X  (t)  ,  (15) 

11       2  2 

where    [X  (t) ,      X  (t) ]   denotes   a   fundamental   system   of 
1       2 

solutions  to  (13) ,    and  C  and  C   are  constants  determined  by 

1       2 

the   initial  conditions  [see,   for   example,   Ref  7].   The 

Y-force  level  has  similar  form.  The  functions  X  (t)  ,      X  (t)  , 

1       2 

Y  (t) t      and   Y  (t)   are  referred   to  by  Taylor  and  Brown  as 
1  2 

General   Lanchester   Functions.    Their    properties   are 

summarized   in   Table   I.   The  functions  X  (t)  and  Y  (t)  are 

1         1 

similar  to  the  hyperbolic  cosine,  while  X  (t)  and  Y  (t)   are 

2         2 

similar   to   the   hyperbolic   sine.   The   General  Lanchester 

Functions  may  be  constructed  by  successive  approximations  or 
by  infinite  series  methods. 
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A.  APPLICATION  TO  POWER  ATTRITION-RATE  COEFFICIENTS 


When  this  solution  method  is  applied  to  the  case  of 
power  attrition-rate  coefficients  (7) ,  the  solution  (15)  may 
be  written  as 

x(t)  =x  [0    (C)  u    (t+C)  -  V    (C)v    (t+C)  ]  - 

0   m,n    m,n         m,n     m,n 

y  Jk  /k   [u    (C)v    (t+C)  -u    v    (C)  ],     (16) 
0   a   b    m,n     m,n         m,n  m,n 

where  u   ,  v   ,  0    and  7    are  referred  to  by  Taylor  as 
tn,  n    m,  n    m,  n      m,  n 

Power   Lanchester   Functions.  They  have  the  properties  shown 

in  Table  II.  For  computational  convenience,   utilization  is 

made  of  the  Lanchester-Clif f ord-Schlaf li  (LCS)  functions 

2k     k"1 
F  (x)  =  Z(x/2)    /  [k  n(j*v)  ], 

v     k=0  j=0 

2k+1      k 
G  (x)  =  £(x/2)      /  [k  n  (j  +  v)  ]. 

v      k=0  j=0 

The  Power  Lanchester  Functions  may  be  expressed  in  terms  of 
the  LCS  functions  through  the  following  relationships: 

u    (t)  =  F  (S(t)) 

m,  n        g 

(m-n)/2 
v    (t)  =  t       G  (S(t)) 

m,n  p 

0    (t)  =  F  (S(t)) 

mrn       p 

(m-n)/2 
V    (t)  =  t        G  (S(t))  r 

mrn  g 


m     n 

where  S  (t)  ~\l  (k  t  )  (k  t  )  t  /  ((m+n+2)/2).    Using   the   LCS 

a     b 

functions,  (16)  may  now  be  written  as 
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x(t)  =  x  fF  (J(O)JF  (J(t))  -  (1+t/C)  (m  Q)/  G  (J(0))G  (J(t)) 

0   p        g  q         p 

(m-n)/2 


-  yja(t=0)/b(t=0)  [  (1+t/C)        F  (J(0))G  (J  (t)  )  - 
0  q        p 

G  (J(0))F  (J(t))  ].  (17) 

p     g 


B.  APPLICATION  TO  OFFSET  LINEAR  COEFFICIENTS 


Similar  techniques  are  utilized  to  apply  this  general 
theory  to  the  case  of  the  offset  linear  attrition-rate 
coefficients  (8)  .  With  these  coefficients,  a  fundamental 
system  of  solutions  to  (13)  is  given  by 

x  (t)=f(t*C),   x  (t)=g(t+C),   y  (t)=F(t+C),  y  <t)=G(t+C) 


where 


2n        n  k  k  4n-k 


f  (t)  =  Z  £  (Vkk  /2)   /  (2n)  !  }  I  B  A  t 

n=0    a  b  k=0n 


2n+1  n  k  k  4n+2-k 


g(t)  =  E  {Wk  k  /2)     /(2n  +  1)!}  EC  A  t 

n  =  0    a  b  k=0n 

2n        n  k  k  4n-k 


F  (t)  =  Z   {  (\/k  k  /2)   /  (2n)  !  }  S  D  A  t 

n  =  0    a  b  k=0n 

2n+1         n+ik  k  4n+2-k 

G  (t)    =     Z   [  (\fk~~ k~  /2)  /(2n  +  1)  !  }   Z  E   A   t 

n  =  0  a    b  k=0n 

k        k  k  k 

The   coefficients    B    ,    C    ,      D    ,      and      E        satisfy      complicated 

n        n  n  n 

recurrence  relations  [see  Taylor  and  Brown  (Ref  13) ].   The 

functions  f  (t)  ,  g  (t)  ,  F  (t)  ,  and  G  (t)  are  referred  to  as 
Offset  Linear  Lanchester  Functions.  As  in  the  case  of  the 
power  attrition-rate  coefficients,  use  is  made  of  the 
following  Auxiliary  Offset  Linear  Lanchester  Functions  for 
computational  convenience: 


18 


oo   2n        n  k  k 
h(p,u)  =  Z  {p   /(2n)  !}  IB  p 

n=0  k=0n 

oo   2n+1  n  k  k 

w(p,u)  =  Z  (p    /(2n+1).  }  z  C  u 

n  =  0  k=0n 

oo   2n         n  k  k 
H(p,u)  =  E  {p  /(2n)!}  Z  D  u 

n  =  0  k=on 

2n+1         n+1k  k 
W(P,u)  =  Z  (p    /(2n+1)!}  Z  E  u  . 

n=0  k=0n 

Utilizing  these  Auxiliary  Lanchester  Functions,  the  X-force 
level  solution  (15)  may  be  written  as 

x(t)  =  x  [H(J(0),j(0))  h(J(t)#j(t))  -  W(J<0)  ,j<0)) 

0 

r(J(t)rj(t))  ] 


-  yjk  /k  £h(J(0)  rj(0))  w(J(t),j(t))  -  w(J(0),j(0)) 

0   a  b 

MJ(t),j(t))  ],  (18) 


where  j  (t)  =  A/(t+C)  and  J(t)  =\/k  k   (t+C)  /2- 

a  b 


The  equations  (17)  and  (18)  are  solutions  to  the 
X-force  level  equation  (13)  for  the  case  of  power  and  offset 
linear  attrition-rate  coefficients,  respectively.  This 
solution  method  of  Taylor  and  Brown  represents  a 
considerable  extension  of  Lanchester  theory  in  the  area  of 
variable  attrition-rate  coefficients.  Previous  analytic 
results  had  been  available  only  in  the  more  limited  case  of 
the  opening  range  of  battle  equal  to  the  minimum  of  the 
maximum  effective  ranges  of  the  two  weapon  systems.  That  is, 
the  initial  fire  effectiveness  of  at  least  one  side  was 
restricted  to  be  zero  [see  for  example  Bef  12]. 

While  the  time  solutions  (17)  and  (18)  appear 
formidable,  their  implementation  is  simplified  if  a  digital 
computer   is   available.    As   noted   earlier  in  this  paper, 
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however,  these  solutions  are  complex  and  involve  new 
mathematical  functions.  This  makes  it  difficult  for  the 
military  analyst  to  perceive  significant  relationships  among 
battle  variables,  and  renders  parametric  analysis  difficult 
and  time-consuming.  This  provides  motivation  for  the 
approximation  technique  to  follow. 
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V.   REDUCTION  TO  LIOOVILLE' S  NORMAL  FORM  AND  APPROXIMATION 

The  argument  ab  (referred  to  by,  among  others,  Taylor 
and  Brown  [Ref  13]  as  the  "intensity  of  combat")  in  the 
constant  coefficient  result  (4)  provides  motivation  for  the 
transformation  £ Ref  14] 


z  =  /\/a(s)b(s)  ds,  (19) 

c0 
where  z  (t=0)  is  denoted  z  .   Applying  this  transformation  to 

0 

(13)  and  (14)  yields 

2    2 
d  x/dz   ♦  {1/2  d/dz  ln£b(t)/a(t)  ]}  dx/dz  -  x  =  0,  (20) 

with  initial  conditions 

1/2 
x(z=z  )  =  x^   and   {[  b  (t)  /a  (t)  ]     dx/dz}    =  -y  , 
0     0  z=z     0 

and 

2    2 
d  y/dz  ♦  (1/2  d/dz  ln[a(t)/b(t)  ]}  dy/dz  -  y  =  0,  (21) 

with  initial  conditions 

V2 
y(z=z  )  =y   and   {[  a  (t)  /b  (t)  ]    dy/dz}    =  -x  . 
0     0  z=z     0 


Taylor  and  Brown  have  shown  that  the  X  force  level 
equation  (13)  may  be  transformed  into  a  linear  second  order 
differential  equation  with  constant  coefficients  if  and  only 
if  (20)  is  a  constant  coefficient  equation.  This  is  the  case 
if  and  only  if 


(1/\/a(t)b  (t))  d/dt  ln[a(t)/b  (t)  ]  =  CONSTANT. 


Equation  (20)  shows  the  significance  of  the   parameters 
a(t)/b(t)   (the   relative   effectiveness   of   the  two  weapon 
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systems)  ,  and  the  intensity  of  combat  \l a  (t)  b  (t)  in 
determining  battle  outcome.  This  is  a  generalization  of  the 
well-known  constant  coefficient  results  to  the  case  of 
temporal  variations  in  fire  effectiveness.  The  significance 
of  these  parameters  may  be  seen  more  explictly  by 
transformation  to  Liouville^  normal  form  (see  Ref  14). 

The  transformation 

1/4 
x(z)  =  X(z)[(a(t)/a  )  /  (b(t)/bQ)]    ,      (22) 

may  be  applied  to  (20)  and  yields  the  normal  form  with  the 
first  derivative  of  the  dependent  variable  removed 

2    2 
d  X/dz   -  [1  +  F(z)  ]X  =  0,  (23) 

with  initial  conditions 


X  (z=z  )  =  x    and   dX/dz(z=z  )  =  -y  \Ja.   /b   -  x  e  , 
00  0       00000 


where 


F(z)  =  P   (z)  /  P(z)      P(z)  =  [H(t)  ]  (24) 


H(t)  =  a(t)  /b(t)      e(t)  =  1/4  \/a(t)  b  (t)  d/dt  InR  (t)  . 

(25) 
2     2 
In  (24),  the  notation  P,,(z)  denotes  d  P/dz  .   Equation  (23) 

is  Liouville's  normal  form  [see  page  23  of  Ref  9]. 

For  the  Y  force  level,  the  transformation 

1/4 
y(z)  ■  I(z)[  (b(t)/b  )  /  (a(t)/a)] 

yields  the  normal  form 

2    2 
d  Y/dz   -  [ 1  +  G(z)  ]Y  =  0,  (26) 

with  initial  conditions 


22 


where 


Y(z=z  )    =   y        and     dY/dz(z=z   )    =  "X Jb  /a     +   ye   , 
00  0  00000 


1/4 
G(z)    =Q«»(z)    /  Q(z)       and      Q  (z)    =  [  R  (t)  ]        =    1/P  (z) 


Expressing  the  transformed  force-level  equation  (23)  in 
the  simple  form 

2     2 
d  X/dz   -  X  =  Q(z)  X, 

variation  of  parameters  may  be  used  to  obtain  the   solution 
to  (23)  as 


X  (z)  =  x  cosh  (z-z  )  -  [y  \j a  /b   +  x  e  ]  sinh  (z-z  ) 

0      zo      00000         0 


♦  /F(j)  sinh(z-j)X  (j)dj.  (27) 


Z0 


This  equation  may  be  expressed  in  terms  of  the  original  time 
variable  t  and  dependent  variable  x  as 

1/4 


x(t)    =   [  (a(t)/a   )/(b(t)/b    ]         {x   cosh  (  /*  (s)  b  (s)  ds) 

o  o  to  o 


-[y  /a/b     +   x   e    ]sinh  (  JVa  (s)  b  (s)  ds) } 
0      0      0  0    0  o 

1/4  1/4 

+  [a(t)/b(t)  ]  d/Js{1/(a(s)b(s))d/ds[b(s)/a(s)  ]        } 

sinh(  /i/a(sj)  b(  j)  d j)  X  (s)  ds.  (28) 

s 
If   the   relative   effectiveness   varies   slowly,   then   the 

relationship   (24)   would   indicate  that  F(z)<<1  so  that  one 

could  drop  the  integral  term  in  (27)  yielding 


X  (z)  =  x  cosh  (z-z  )  -  [y  va/b   +  x  e  ]  sinh  (z-z  ), 
0         0       0   0   0     0  0  0( 

which   is   the  Liouville-Green   approximation  [Ref  11].  The 

approximation  may  be  expressed   in   terms   of  the   original 
variables  t  and  x  as 
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x(t)    «  [  (a(t)/a  )/b(t)/bQ] 


1/4 


(x   cosh  (  /Va  (s)  b  (s)  ds) 
t°  0 


-  [y   va/b     ♦   x   e    lsinh  (  /ya  (s)  b  (s)  ds) } 
0      0      0  0    0  o 


(30) 
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71.    APPLICATION  TO  POWER  AND  OFFSET  LINEAR  COEFFICIENTS 

The  Liouville-Green  approximation  may  now  be  applied  to 

the    constant-speed   attack    model  with  attrition-rate 

coefficients  (7).  For  the  reader's  convenience,  these 
coefficients  are 

m  n 

a(t)  =  Jc  (t+C)     and   b  (t)  =?  k  (t+C)  .      (7) 

a  b 

This  model  assumes  that  the  maximum  effective  ranges  of  the 
two  opposing  weapon  systems  are  equal  (no  offset) .  The 
opening  range  of  battle  is  taken  to  be  within  this  maximum 
effective  range  (i.e.,  C>0).  Applying  the  Liouville-Green 
approximation  (29)  ,  the  time  solution  to  the  X  force  level 
equation  in  terms  of  the  original  time  variable  t  and 
dependent  variable  X  becomes 

(m-n)/4        5 


x(t)  =  (1+t/c)        (x  cosh./Wa  (s)b(s)  ds 

0     0 

, (m-n)/2  -(m+n+2)/2 

-  [y  Vk  /k  C       +  (x  (m-n)/4  k  k  )C  ] 

0   a   b         t  0         a  b 


sinhjVa  (s)b  (s)  ds} ,  (31) 

0 


where 


f\la.  (s)b(s)  ds  =  (Vkk  /(m+n+2)/2) 
0  a  b 

(m  +  n  +  2)/2   (m+n+2)/2 
£  (t+C)         -C         }  . 


For  future  reference,  equation   (24)   becomes   in   this 
case 

2  2 
F(z)  =  [  (m-n)  (3m+n+4)  ]  /  [4(m+n  +  2)  z  ],      (32) 

or  in  terms  of   the   original  time   variable   t   and   range 
parameters  (11) : 
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ffl  n 
F(t)  =  [B  R  (m-n)  (3m+n+4)  ]  / 
a  b 

nn-n  m+n+2 

[16a  b  (vel)    <t+(R  -R  )/vel)      ],       (33) 
0  0  a   0  ' 

where  vel  =  constant  velocity  of  attack. 

We  can  also  apply  the  Liouville-Green  approximation  to 
the  case  of  the  offset  linear  attrition-rate  coefficients 

a(t)  =  k  (t+C)    and   b  (t)  =  k  (t+C+A)  .      (8) 
a  b 

In      this      case,      the      time      solution      to      the      X   force   level 
equation   becomes 

1/4 


x(t)    =   [  (1+(A/C)/(1+(A/(t+C)))  ]  '      {x    cosh/\/a(s)b(s)  ds 

0  0 

-   [  (jjk   /k J/(1+(A/C))     ♦     (x     (A/C))/ 

0      a     b  0 


3/2 


(4  k   k   C    (1+(A/C))         )  ]sinh/Va(s)b(s)  ds}  ,  (34) 

a   b  0 


where 


JVa(s)b(s)  ds  =    (A      k   k  /8)  {  (1+ (2  (t  +  C)/A)  )     (1+ (2  (t+C) /A) )    -1 

0  a   b 

2 
-    ln[  1+ (2  (t*C) /A)     +       (1+ (2  (t+C) /A))    -1]} 


For   future   reference,    equation       (24)       becomes      in      this 
case 

F(t)     =    R   fi    (R    -R    )[12(t+(R   -R    ) /vel)     +    7  { (R    -R    )/vel)  ]   / 

aBBa  aO  Ba 

3  3  3 

{16a  b   vel    (t+ (R   -R    )/vel)     (t+ (R    -R   )/vel)     }.       (35) 
0   0  a      0  b      0 
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VII.   DISCUSSION 

A  computer  program  was  developed  in  FORTRAN  to  produce 
plots  of  the  X  force  level  versus  range  between  opposing 
forces  for  the  constant-speed  attack  model  described 
earlier.  Both  the  power  and  offset  linear  attrition-rate 
coefficient  cases  (equations  (31)  and  (34))  were  considered. 
These  approximate  solutions  were  compared  to  the  exact 
results  given  by  Taylor  and  Brown  [Ref  13]  to  evaluate  the 
accuracy  of  the  Liouville-Green  approximation  to  equations 
with  these  coefficients. 

A.  POWER  ATTRITION-RATE  COEFFICIENTS 

For  the  case  of  power  attrition-rate  coefficients,   the 

maximum   effective  ranges  of  the  two  weapon  systems  were  set 

equal  to  2000  meters,  with  an  opening   range   of   battle  of 

1250  meters.  The  fire  effectiveness  of  the  two  sides  at  zero 

range  separation  was  a  =.06   X  casualties/min-Y   unit,   and 

0 

b  =.6   Y   casualties/min-X  unit.   The  velocity  of  attack  was 
0 

constant  at  five  miles   per   hour,   and   the   initial   force 

levels   were   X  =   10   and   Y  =  30.   These  parameter  values 

correspond  to  those  used  by  Taylor  and  Brown  in  their   paper 

[Ref   13].    The   exponents   m   and   n   in   (7)   were  varied 

throughout  the  range  explored  in  Ref   13.   The   results   are 

shown   in   Figures   3  through  10,  where  each  figure  contains 

the  approximate  time  solution  and   the   corresponding   exact 

analytical  solution  obtained  by  Taylor  and  Brown.   As  even  a 

cursory  examination  reveals,  the  approximation   appears  to 

conform   to   the   established  results  only  in  certain  cases. 

Some  possible  explanations  for  this  observed  deviation   will 

now  be  considered. 

Recall  that  the  integral  term   in   equation   (27)   was 
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dropped  (assumed  to  be  approximately  zero)  in  making  the 
Liouville-Green  approximation.  The  validity  of  this 
assumption  appears  to  be  linked  strongly  to  the  form  of 
F  (z) ,  which  for  the  present  case  is 

0  n 
F(t)  =  [R  E  (m-n)  (3m+n+4)  ]  / 
a  b 

m+n  m+n+2 

[16a  b  (vel)    (t+  (^"Rq)  /vel)      J.       (33) 

Hhen  the  exponents   m  and  n   are  equal,   F(t)=0/  and  by 

considering   equation    (27) ,    one  can    see   that   the 

approximation  is,   in  fact,   the  exact  solution.   This  is 

verified  trivially   by  Figures   3  and  4,   for  which  the 

exponents  are  equal.   It   should  be   noted  that  in  this 

situation,  the  model  reduces  to  the  case  of  a  constant  ratio 
of  attrition-rate  coefficients. 

When  the  exponent  m  is  greater  than  n,  F  (t)  will  always 
be  greater  than  zero,  and  from  (27)  it  follows  that  the 
exact  solution  is  greater  than  the  approximation.  This 
situation  is  illustrated  by  Figures  6  and  10.  The  converse 
also  holds.  When  m  is  less  than  n,  F (t)  is  always  less  than 
zero,  and  the  approximation  is  greater  than  the  exact 
solution.  This  is  shown  in  Figures  5  and  9. 


Again  examining  (33)  ,  the  parameter  la  b  I ,  which  is  the 

1  0  0 

product  of  the  opposing  forces'  fire  effectivenesses  at  r=0, 

appears   to   be  significant.   is  la  b    increases,   F(t) 

0  0  ' 

decreases,  making  the  integral  term  in  (27)  smaller,  so   one 

would  expect  the  approximation  to  be  "closer"  to  the  exact 
solution.  This  did  occur  in  the  cases  considered  (see  for 
example  Figure  7) . 


B.  OFFSET  LINEAR  COEFFICIENTS 
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For  the  case  of  the   offset   linear  attrition-rate 

coefficients  (8)  ,  the  maximum  effective  range  of  the  Y  force 

weapon   system   (R  )   was  set  equal  to  1500  meters,  with  the 

a 

rest  of  the  parameter  values  equal  to  those  in   the   power 

coefficient  case.  The  maximum  effective  range  of  the  X  force 

weapon  system   (H  )   was   varied   throughout   the   range 

b 

considered   in  Ref  13,   The  results  are  shown  in  Figures  11 

through  15  where,  as  before,  each  figure  compares  the 
corresponding  solutions,.  As  was  the  case  with  the  power 
coefficients,  the  solutions  are  seen  to  conform  only  in 
certain  cases. 


For  the  linear  coefficients,  it   was   shown   that   F  (z) 
could  be  written  in  terms  of  the  time  variable  t  as 

F(t)  =  R  R  (R  -R  )[12(t+(H  -R  ) /vel)  +  7  (  (R  -R  )/vel)  ]  / 

aBBa         aO  Ba 

3  3  3     ,   % 

{16a  b  vel  (t+ (R  -R  )  /vel)  (t+ (R  -B  )/vel)  }.    (35) 
0  0         a   0  b   0 

When  the  aaximum  effective  ranges   of   the   opposing   weapon 

systems   (R   and   R  )  are  equal,  F (t)  =  0,  and  the  integral 
a        b 

term  in  (27)  is  zero.   This  again  reduces   to  the   constant 

attrition-rate-ratio  case  as  can  be  seen  from  (8)  with  A=0. 
Thus  one  would  again  expect  agreement  of  the  exact  and 
approximate  solutions.  This  was  found  to  be  true  as  shown 
by  Figure  11,  where  both  maximum  effective  ranges  equal  1500 
meters. 


As   was   the   case   for   the   power   coefficients,   the 

parameter  la  b    appears   in  the  denominator  of  F(t),  hence 
'0  0' 


the  two  solutions  should  be  closer  as  a  b 

1  0  0 
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increases.   This 


was  found  to  be  the  case,  as  illustrated  by  Figure  15. 


C.   ANALYTICAL  CONSIDERATIONS 

The  examinations   of   both  the  power  and  offset  linear 

attrition-rate  coefficient   cases   have  centered   on   the 

parameter  F  (z)  in  the  solution  (27) .  While  F  (z)  does  appear 

to  be  the  driving  factor  in  the  integral  term  in  (27)   (thus 

in   the   validity   of   the  approximation) ,  it  should  be  made 

clear  that  the  other  factors  in  the  integral  term  can  not  be 

ignored.    While   a  mathematical  analysis (for  example  of  the 

hyperbolic  sine  term)  would  be  formidible,  examples   can  be 

found   which   illustrate  that   F (z)   can  not  always  be 

considered  alone.  One  such  example  is  shown  in  Figure  8.  The 

fire   effectiveness   of    the   two  forces  at   r=0   were 

interchanged  (i.e.,  a  =  .6  and  b  =.06  instead  of  vice   versa) 

0         0 

in   the   power  attrition-rate   coefficient   case.   From  the 

appropriate  form  of  F  (z)  ,  equation  (33),  one  can  see  that 
this  change  has  no  effect  on  the  numerical  value  of  F (t) 
throughout  the  time  range.  let  the  significant  difference  in 
results  can  be  seen  in  Figure  8  when  compared  to  Figure  6. 
Currently,  there  is  no  convenient  methodology  to  analyze  the 
implications  of  such  discrepancies. 

Toward  this  end,  theoretical  analysis  in  the  area  of 
Olver's  "Error  Bounds  for  the  Liouville-Green  Approximation" 
[ Ref  11]  could  be  very  beneficial.  While  they  sound 
promising,  Olver1 s  results  can  not  be  directly  applied  to 
variable-coefficient  Lanchester-type  differential  eguations. 
If  available  in  a  convenient  form,  such  error  bounds  could 
be  evaluated  for  specific  situations  to  investigate 
tolerable  error  limits  prior  to  invoking  the  approximation. 

From  the  above  analytic   considerations,   one   can   see 
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that  while  one  may  be  able  to  explain  some  of  the  error 
incurred,  this  is  small  comfort  to  the  analyst  attempting  to 
make  practical  use  of  the  Liouville-Green  approximation.  It 
has  been  shown  that  the  approximation  appears  to  yield 
credible  results  only  in  certain  limited  circumstances. 
Furthermore,  methodology  is  not  presently  available  to 
predict  the  error  which  will  be  incurred  in  a  specific 
situation,  except  for  the  trivial  case  in  which  the  integral 
term  in  (27)  goes  to  zero. 
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VIII.    SOHHABY 

This  thesis  has  examined  the  "adequacy"  of  the 
so-called  Liouville-Green  approximation  to  the  solution  of 
variable-coefficient  Lanchester-type  equations  of  modern 
warfare.  Although  this  approximation  only  involves 
"elementary"  functions  and  has  an  intuitively  appealing 
form,  it  was  unfortunately  found  that  this 
Liouville-Green-Lanchester  approximation  is  not  consistently 
reliable  for  ^estimating  force  levels.  Specifically,  the 
approximation  was  applied  to  two  forms  of  attrition-rate 
coefficients:  (I)  power  attrition-rate  coefficients 
(modelling,  for  example,  the  same  maximum  effective  ranges 
for  the  opposing  weapon  systems) ,  and  (II)  offset  linear 
attrition-rate  coefficients  (modelling  different  maximum 
effective  ranges) .  It  was  seen  that  the  approximation  was 
made  by  dropping  the  integral  term  in  equation  (27) ,  hence 
exact  results  are  achieved  only  when  this  term  equals  zero. 
This  situation  was  found  to  occur  only  when  the  ratio  of 
attrition-rate  coefficients  was  constant.  For  the  power 
attrition-^rate  coefficients,  this  implies  that  the  exponents 
in  (7)  are  equal,  while  for  the  offset  linear  coefficients 
this  implies  the  offset  parameter  A  in  (8)  is  zero  (which 
occurs  when  the  maximum  effective  ranges  are  equal) . 
Analytic  considerations  showed  that  while  certain  initial 
parameter  combinations  produced  credible  results  from  the 
approximation,  there  was  no  methodology  available  to  predict 
the  error  incurred  in  the  general  case.  It  was  found  that, 
as  a  general  rule,  the  further  one  moves  from  the  constant 
coefficient  ratio  case,  the  larger  the  error  incurred  by 
invoking  the  Liouville-Green  approximation. 
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TABLE   I. 


Properties  of   the  General  Lanchester  Functions     x. ,   x.,   y   ,   y~, 


1.  dx^dt  =  \/\  a(t)y2 


2.  dx2/dt  =   \/\  a(t)yi 


3.  dyx/dt  =  /ka7k^b(t)x2 


4.  dy2/dt  =  ^J\  b(t)Xl 


5.  Xl(t)yi(t)    -  x2(t)y2(t)   =1  Vt 
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Table  II.   Properties  of  the  Power  Lanchester  Functions 

u   ,   v   ,   U   ,  and  V 

m,n    ra,n    m,n         m,n 


i.  dum  /dt  =  or  tk  tm]v    <t) 

m,n  b     a       a  ra,n 


2.                       dv       /dt  =  /k,/k  [ktlU        (t) 

m,n  b     a  a          tn,n 

3«                                    dUm    n/dt    =  &JK  ^K^r*    Jl) 

m,n  a      b  b          m,n 

4.                      dV       /dt  =  /k  /k,  [k,  tn]u       (t) 

m,n  a     b  b          m,n 


5.  u        (t)U        (t)   -  v        (t)V        (t)    =1  Vt 

m,a  m,n  m,n  rn,n 


6.  u        (t=0)    =  U        (t=0)    =  1 

m,n  m,n 


7.  v        (t=0)    =  V        (t=0)    =  0 

m,n  m,n 


m+1 


8.  u     m(t)    =  Um      (t)    =  cosh    (/k~k7  t      V(nH-l)) 

m,m  m,m  a  b 


m+1 


9.  v        (t)   -  V        (t)    =  sinh   (/k~k7  t     V(nrJ-l)) 

n,m  m,m  a  b 


Properties   6  and  7  only  hold  for     m,n    >  -1. 
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Figure  3.  Power  coefficients  with  m=n=l. 

For  Figures  3  through  10,  R  =R[=2000m .  , 
Xq=10,  y0=30,  aQ=.06,  and  bn=.b  unless 
otherwise  specified. 
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Figure  4.  Power  coefficients  with  m=n=2 
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Figure  5.  Power  coefficients  with  m=2,  n=3 
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Figure  6.  Power  coefficients  with  m=3,  n=2 
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Figure  7.  Power  coefficients  with  m=3,  n-2,  a  =B  =.6 
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Figure  8.  Power  coefficients  with  m=3,  n  =  2,  a  =  .6,  B  =  .06 
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Figure  9.  Power  coefficients  with  m  =  2,  n=2.1 
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Figure  10.  Power  coefficients  with  m=2.1,  n  =  2 
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Figure  11.  Linear  coefficients  with  R  =  1500m. 

For  Figures  11  through  16, B  Ra=1500m., 
x0=10,  ^0=30'  a0=-06>  and  bQ=.6  unless 
otherwise  specified. 
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Figure  12.  Linear  coefficients  with  R  =  1600m 
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Figure  13.  Linear  coefficients  with  R  =  2000m 
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Figure  14.  Linear  coefficients  with  R  =  2500m 

B 


48 


a 
a 


o 

N 


O 
Ci 


a 

\u  a 

Q 
U. 


a 

o 


Exact 


Approx . 


2.50 


5.00 
RANGE 


1X10     ) 


7.50 
2 


I 

10.00 


12.50 


Figure  15.  Linear  coefficients  with  R  =  3000m 
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Figure  16.  Linear  coefficients  with  R  =2500m.  ,  a  =b  =.6 
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